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Abstract 

X-ray four-wave mixing signals generated in the ki + k2 — phase-matching direction are 
simulated for Nls transitions in para-nitroanline and two-ring hydrocarbons disubstituted with 
an amine and a nitroso groups. The two-dimensional x-ray correlation spectra (2DXCS) provide 
background-free probes of couplings between core-electron transitions even for multiple core shells 
of the same type. Features attributed to couplings between spatially-separated core transitions 
connected by delocalized valence excitations provide information about molecular geometry and 
electronic structure unavailable from linear near-edge x-ray absorption (XANES). 
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I. INTRODUCTION 



Near-edge x-ray absorption spectroscopy (XANES) provides a powerful frequency-domain 
probe for electronic structure of molecules. 11] Transitions from the ground state to bound 
core-excited states appear as resonances in the absorption spectrum below the ionization 
edge. Because of the compactness of core shells, the positions and intensities of XANES 
peaks arising from a given shell carry information about the electronic structure in its vicin- 
ity. XANES carries characteristic signatures of the electronic environment of the absorbing 
atom. If two atoms are spatially well separated, their contribution to XANES is essen- 
tially additive. The molecular structure can then be elucidated by identifying the signatures 
of various functional groups in the total XANES spectra. This additivity, known as the 
building-block principle of XANES makes it insensitive to electronic-structure variations 
away from the absorbing atoms as well as to subtle differences in molecular geometry. 

Coherent nonlinear x-ray spectroscopies can overcome these limitations and extend 
the XANES capabilities towards more detailed probes of electronic and molecular struc- 
ture. Coherent nonlinear measurements performed with infrared and visible femtosecond 
phase-locked pulse sequences, can enhance desired spectral features, eliminate certain line- 
broadening mechanisms, and detect interferences between specific quantum pathways con- 
tributing to the optical response. The ongoing development of high- 
harmonic generation (HHG) and forth-generation synchratron sources based on the x-ray 



free-electron laser (XFEL) (see Refs. 
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and references therein) provide first steps 



towards the realization of coherent nonlinear measurements in the x-ray domain. These will 
require multiple x-ray pulses with controlled timing, phases and sufficient intensity. 

As these sources continue to develop, one may rely on theoretical simulations to design 



and evaluate possible nonlinear techniques. Earlier studies focused on u 



sorption and scattering in systems prepared by an optical pulse 
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and Mukamel studied frequency- domain all-x-ray four- wave mixing 18|, ll9|] . Pump-probe 



trafast x-ray ab- 



171]. Tanaka 



is the simplest time-domain nonlinear experiment. This incoherent technique requires two 



pulses with variable time 
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delay but no contro 
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2l| or infrared 



22 



over the phases. Combinations of optical 
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24| ) and x-ray probe as well as x-ray 



pump (either visib 
pump/x-ray probe 

four-wave-mixing techniques which require up to four x-ray beams. 2m, l28|] These offer a 



261 ] have been studied. We consider on attosecond phase-coherent 
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much higher degree of control of the observed dynamical processes and could result in quali- 
tatively new information unavailable from any other technique. In an earlier study, we have 
examined the ki = —ki + k2 + k-^ 2DXCS signal of aminophenol obtained by varying two 
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29l The simulated 



delay periods in the coherent x-ray four- wave mixing measurement, 
two-color 2DXCS signal where two pulses are tuned to the N K-edge and the other two to the 
O K-edge was shown to be highly sensitive to the coupling of the spatially- and spectrally- 
separated core transitions. Distinct off-diagonal cross peaks appear due to the interference 
among quantum pathways that involve only singly core-excited states (excited-state stimu- 
lated emission [ESE] and ground-state bleaching [GSB]) and pathways that involve singly 
and doubly core-excited states (excited-state absorption, ESA). If the frequency of a given 
core-shell transition is independent of whether another core-shell is excited, the ESA con- 
tribution interfers destructively with the ESE and GSB and the cross peaks vanish. The 
coupling between two transitions results in a distinct 2DXCS cross-peak pattern. In con- 
strast, XANES of two independent transitions is exactly the sum of the individual transition. 
Since the coupling between two core transitions depends on the distance between the two 
core shells as well as the electronic structure in their vicinity, 2DXCS cross peaks carry a 
wealth of qualitatively new information beyond XANES. 

The simulated kj signal of aminophenol has a simple structure because the 100 eV separa- 
tion between the N Is and O Is transitions is much larger than the assumed pulse bandwidths 
(< 10 eV). The two-color 2DXCS signal thus involves transitions from both cores and the 
resulting spectrum contains no diagonal peaks. Here we focus on a homonuclear 2DXCS 
signal in systems with multiple core shells of the same type. In this case, both transitions in- 
volving two different or the same core-shells contribute to the signal since the chemical shifts 
(a few eV) are smaller than the pulse bandwidths and the 2DXCS diagonal and cross peaks 
spectrally overlap. Due to interference, the latter are usually weaker. A higher spectral res- 
olution is thus required to separate the cross peaks and extract the couplings. The kj signal 
of nucleobases and their pairs was shown to be dominated by a strong GSB contribution 
arising from transition of imine N Is into the vr orbitals of the heterocycle. 30|] 

In this paper, we show that the single- resonance contributions can be eliminated by 
monitoring the 2DXCS signal in the kjjj = ki + fc^ — ^3 direction. This technique 31 1 
corresponding to double-quantum coherence in NM R|32 | w as already predicted to show high 
sensitivity to exciton coupling in the infrared 3^, 34, 35| and the visible {sgI, 



37|. Within 
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the rotating-wave approximation, only two ESA pathways contribute to this signal, both 
involving doubly core-excited states. When two core transitions are independent the two 
pathways intefere destructively. This signal thus contains only features induced by the 
coupling between core transitions. 



In Section [IT] we employ the response-function formalism 38|] to derive the sum-over- 
states expression for the km signal. The kj and kni signals of a model four-level system 
with and without coupling between the two core transitions are compared in Section lllli In 
Section IIVI we present the Nls XANES and km 2DXCS signals of benzene, stilbene, and 
biphenyl disubstituted with the amine and nitroso groups (Fig. [1]). The relevant core-excited 
states are described using singly- and doubly-substituted Kohn-Sham determinants within 
the equivalent- core approximation |29:]. The results are summarized in Section IVl 

II. SUM-OVER-STATES EXPRESSIONS FOR THE fci + fca - ^3 SIGNAL 

The most general four-wave-mixing experiment uses a sequence of four soft x-ray pulses 
(Fig. EI). Possible experiments with fewer pulses are discussed below. 

4 

E{r, t) = J2 hS,{t - r,)e^^^™^(*--^) + c.c] (1) 
j 

Here, Sjir) is the complex (positive frequency) temporal envelope of the j'th pulse, ej is 
the polarization, kj is the wavevector, and uj is the carrier frequency. Since the size of 
the N core shell (~ 0.03 nm) is much smaller than the N K-edge wavelength (~ 3 nm), we 
can safely use the dipole approximation to describe the interaction of the x-ray field with a 
molecule, 

4 

Hint = -Yl [f'^^'J^^ - r,)e^'=^^-*-^(*--^) + C.C.] . (2) 

Here, fi^ = (e^ • /i)e**'j^^ and Vj is the laboratory-frame position of the core-shell interacting 
with pulse j. The heterodyne kjjj signal is recorded as a function of the three delays ^1,^2, h 
between consecutive pulses (Fig. [2]) 

//•oo 
dr / dtP(3)(r,t)^:(t)e''=*^-^-** (3) 

where the induced polarization P^^^ is calculated via third-order time-dependent perturba- 
tion expansion in the field [27|]. The spectrum will be displayed as the Fourier transform of 



Siii with respect to ts and ^2, holding ti fixed 







We consider delays {1 is >tj > 10 fs) longer that the pulse durations (Tj < 1 fs) so that the 



39|with Ai = A2 = +1, 



system interacts with each pulse sequentially. Using the result of Ref. 
A3 = —1, the signal is given by 

C /-n n + ^ ^i;^ "^4 (^2 + UJi - + Ua.^^^ - UJa2b2 - ^axbx)£z (^aafe ~ ^^s) 

bni[}h,^h,ti) = -77^ — ^ 

X ~ 



(^2 — ^a2b2 — ^aibi + ^1^0262 + aib^) 
X ([[[A''Aa363]'Aa262],AU]>' (5) 

where a and E'a are, respectively, the molecular eigenstates and their energies. Here, 
ujab = Ea — Eb is the frequency and Tab is dephasing rate of the transition between states a 
and b; flab = \^a)l^ab{^b\ and S+{uj) = [<?7(cu)]* = J dr£,(r)e^-^ 

The spectral bandwidth of S^i^^ — ujj) is limited to [ct; — cUj] < l/Tj. Only terms where 
i^aibi corresponds to a transition from the ground to a core excited and ujb2a2 to a transition 
from singly- to doubly core excited states thus contribute to the signal. Two terms (Fig. [3]) 
then contribute to Eq. (I5j), 



Slll{^3, ^2, t 

= E 



_ ^si^ge + UJ3)S^{uJef + UJi)£t{^fe' " ^2)gl^(^e'g - UJj) fJ,%fJ,tgfI%, fij, f _i(^^,^_^^)4^_r^,^ 

(6) 

For simplifying the description we have considered an ideal four pulse experiment. Eq. ([6]) 
shows explicitly the roles of the various control parameters. The pulse envelopes £j select 
the core/valence excitations allowed within their bandwidths. The ^2 and ^3 resonances 
show the core excitations during ^2 and t^. Core-exciton dephasing takes place during ti. 

Fewer pulses may be used in practice. The first two pulses may be the same, setting ti = 0. 
By frequency dispersing the signal, beam ^4 can be a long cw pulse and the information 
about gathered in the frequency domain. Thus the experiment may be carried out using 
two ultrashort pulses. 



III. 2DXCS -ki + k2 + k-i AND fci + ^2 - k^. SIGNALS OF MODEL COUPLED AND 
UNCOUPLED FOUR-LEVEL SYSTEMS 

In order to demonstrate the sensitivity of 2DXCS signals to the couphng between transi- 
tions, we have simulated the one-color {ujj = uj) kj and km signals for two model systems 
of uncoupled (Fig. H]) and weakly-coupled (Fig. [5]) core transitions. We assumed broad 
bandwidth S'^^Ufe^ — cu) ~ £^{ujeig — u) ^ 1 for simplicity. 

Using Eqs. (10)-(12) of Ref. |29i, the one-color —ki + ^2 + ^3 signal is given by 

Sj{n3, 0, Qi) = sf'^'ins, 0, fii) + ^^(^3, 0, Qi) + sf'^iQs, 0, (7) 



= \l^gei I \f'gei \ [" 



— Uej^g + iTe-^^g Q-^ — UJe2g + i^e2g'^ ^ ^ 



{Ql + UJe^g + iTe^g){yL^ - Ue^g " A iF/eJ 

/^ge2/^e2/(/^/ei/^eig "I" /^/e2/^e2g) 

(fil + Ue^g + iFe29)(J^3 " ^e^g - ^ + fes) 

where A is the anharmonicity of the single to double transition frequency 



(9) 



^fg = ^eig + i^e2g + ^ (lO) 

The kj 2DXCS of a four-level system in general consists of two diagonal peaks at {—^i = 
uJej^g,Q3 = Ueig) and {—0,1 = uje2g,03 = ^e2g) ^"^^ off-dlagoual cross peaks at (— fii = 
(x^eig, ^3 = ^^625) and (— i^i = ^^625) ^3 = ^e^g)- If two transitions are decoupled (Fig.HJ middle 
column), the spectrum is additive and the ESA contribution to the cross peaks cancels the 
ESE and GSB contributions and the cross peaks vanish. However, since the ESA term 
does not contribute to the diagonal peaks, they remain finite. If the two transitions are 
weakly-coupled (Fig. middle column), the ESA contribution is red shifted by A resulting 
in two-lobe cross peak line shape. However, due to the destructive interference between 
the ESE/GSB and ESA terms, the cross peaks are much weaker than the diagonal peaks. 
Thus, for both decoupled and coupled transitions, high spectral resolution is necessary to 



distinguish the cross peak and extract the information about the couphngs from the kj 
spectrum. 

Using Eq. ([H]), the one-color fci + ^2 — ^3 2DXCS signal in the broadband limit is recast 

as 

n2-UJfg + iTfg 

^3 ~ ^eig + "i^eig ^3 ~ ^e2g + '>'^e2g 

A^geiA'-ei/ f^ge2f^e2f "j /-^-^\ 

^3 - tUe23 - A + iF/ei ^3 - ^^e^g - ^ + fe2 

The spectrum of our four- level system thus consists of two pairs of peaks at {Q3 = uje^g-, ^2 = 
ujfg) and (^3 = uje2gi^2 = ^fg)- Both contributing pathways involve doubly-excited states, 
and the spectrum thus provides single-resonance-free probe of the couplings between the 
transitions. If the two core transitions are decoupled (Fig. HI right column), then ESA2 = 
- ESAi and the total signal vanishes. If the two core transions are weakly coupled (Fig. [5l 
right column), the ESA2 contribution is red shifted by A, and the total spectrum exhibits 
two peaks each having the two-lobe structure similar to the ki cross peaks. The splitting 
between the lobes is approximately equal to the anharmonicity. 



IV. NUMERICAL SIMULATIONS 



2DXCS signals in molecules depend on states with two core-electrons excited. In all cal- 
culations, we used the singly and doubly-substituted Konh-Sham (KS) determinants in the 
equivalent- core approximation (refered to as DFT/ECA) to describe the necessary singly 
and doubly core-excited states. The expressions for the transition frequencies and dipole 



moments within this approximation were presented in Ref. |29|. Doubly-substituted determi 



nants are necessary to describe states in which the two core electrons are promoted to orbitals 
higher than the HOMO. The KS orbitals for the original and equivalent-core molecules where 
obtained with the combination of Becke's exchange 401] and Perdew's correlation 411] func- 
tionals and a combined basis set of Gaussian-type atomic orbitals, whereby an extensive 
IGLO-III 42] set was used on N and a moderate 6-311G** setj43| was used on all other 



atoms. The orbitals and their energies were calculated with the Gaussian 03 electronic 
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structure code. 



This computational protocol was tested by comparing the simulated Nls XANES of 
para-nitroaniline to experimental Nls inner-shell electron e nerg y loss spectroscopy (ISEELS) 



45 



The experimental ISEELS 



of aniline, nitrobenzene, and para-nitroaniline from Ref. 
are displayed on Fig. (left panel). Under the experimental conditions in Ref. |45| (high 
impact energy and small scattering angle) the ISEELS is expected to resemble XANES. 
The exprimental pre-edge ISEELS of para-nitroaniline (Fig. [6], right panel) consists of two 
weaker amine peaks at 401.6 eV and 402.6 eV, a strong nitroso peak at 403.8 eV, and amine 
peak at 405.0 eV. The measured core ionization potential is 406.0 eV. 

The absorption edge {i.e., the frequency of the lowest transition) is given within the 
DFT/ECA as the energy difference between the core orbital in the original molecule and the 
HOMO energy in the equivalent-core molecules. This estimate however neglects the effect of 
the core-shell ionization on the remaining core electons. The relaxation among core electrons 
does not significantly affect the valence electrons, hence, its effect on XANES is limited to 
a shift of the entire spectrum. The core relaxation as well as relativistic effects contributing 
to the core-transition frequency can be corrected for by comparing with experiment. The 
nitroso Nls absorption edge in para-nitroaniline calculated within the EGA underestimates 
the experimental one by 21.7 eV. We further assumed that the core relaxation effects have 
the same magnitude for all system considered as well as the amine N Is edge and thus shifted 
all the EGA core-transition frequencies by this value. 

The main features of NA Nls XANES are qualitatively reproduced by the DFT/EGA 
method. The calculated XANES of para-nitroaniline (NA) overestimates the splitting be- 
tween the lowest amine and the nitroso transitions (marked A and G in Fig. [6]) predicting 
it to be 3.1 eV compared to experiment (2.2 eV). It also overestimates the second amine 
peak (marked B in Fig. [6]) intensity. Also the experimental amine peak at 405.0 eV is not 
reproduced. Instead, the calculated XANES features weak peaks at 404 eV and 406 eV. In 
the the experimental spectrum these peaks may be covered by the strong nitroso peak. 

The calculated XANES of 4-nitro-4/-aminestilbene (NASB, Fig. El right panel) and 4- 
nitro-4/-aminediphenyl (NADP, Fig. [6l right panel) are similar to NA and consist of an 
amine peaks at around 401.0 eV, stronger amine peaks at around 402.6 eV and the strongest 
nitroso peaks at arount 404.0 eV. 

Figure [7] shows the relevant EGA orbitals that give rise to the described features in 
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the calculated XANES. The first amine peak in all three molecules involves promoting the 
Is electron to the vr* orbital of the conjugated system. Note that within the EGA, the 
core-hole potential is described by increasing the nuclear charge by 1, keeping the core- 
shell doubly-occupied, while an extra valence electron is added to describe the promoted 
core electron. The lowest EGA orbital is thus occupied and its shape resembles a bonding 
vr orbital rather then antibonding one. Two-ring NADP and NASB have another '7r*-like 
orbital that gives rise to the very weak peak at 402.0 eV in the XANES spectra. The 
second strong amine peak arises due to the excitation of the amine Is electron to the a* 
orbital of the amine group. Finally, the strongest XANES feature at 404.0 eV arises due to 
transition of the nitroso Is electron to the vr* orbital, which mostly localized on the NO2 
group. The intensity of the 402.5 and 404.0 peaks can thus be explained by the fact that 
the corresponding EGA orbitals are strongly localized on the respective group. 

We next turn to the Nls km signal obtained with four pulses of the same frequency 
ujj = 402. 6e\^, the same temporal envelope, and linear polarization ej = e^. This one-color 
signal of a sample of randomly oriented molecules is given at ti = by 

Sillies, f22, tl = 0) = ^ £~{(^ge + Uj)£~{uJef + Uj)£^{uJfe' - Uj)£^{uJe'g - Uj) 

e,e'J 

where (/i/e'Aie'g/ige/We/) is calculated as described in Appendix Rl 

We assume 100 attosecond rectangular pulse envelopes £j{uj) with 6 eV bandwidth around 
the carrier frequencies {i.e., £j(uj—ujj) = 1 for [ct;— cuj] < 3 eV and zero otherwise). The same 
dephasing rate Tab = 0.3 eV is assumed for all transitions. Transitions to the continuum lie 
outside the chosen bandwidth and are neglected in the simulations. 

km 2DXGS of NA, NADP and NASB are shown in Fig. [81 There are three core-excited 
states with significant dipole strength in the XANES of NA: two due to excitation of the 
amine core electron (states A and B) and one due to the excitation of the nitroso core electron 
(state G). The corresponding features in NA (Fig. [8], left panel) can be identified by their 
position along the ^2 axis. The strongest feature is attributed to the doubly excited state 
corresponding to states B and G. The strong dipole coupling of each state results in a very 
strong individual ESAi and ESA2 contributions. The negative inteference between the two 
pathways leads to the two peaks. The weaker dipole strength of the transition from to ground 



'3 ^ ^en -r i^eq 



1 

1^3 - iOfe + iTfe 



[12) 
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state to state A results in a weaker individual ESAi and ESA2 contributions. However, due 
to the coupling between the A and C transitions, the transition dipole moment changes 
sign and the ESAi and ESA2 pathways intefere constructively, resulting in a strong overall 
signal. Comparison of the equivalent-core orbitals describing the promoted core electron 
in the singly and doubly excited states (Fig. [7] and [HI left panels) explains the difference 
between the two sets of resonances. In state A, the promoted core electron is delocalized, 
hence states A and C states are strongly coupled, and the dipole coupling between the A 
states and the doubly excited state has opposite sign to the coupling between the C state 
and the doubly excited state. In contrast, state B is localized at the amine group, hence, 
the coupling between states B and C is weaker, and the overall signal is much smaller than 
each of the individual contributions. 

Similar to NA, there are three core-excited states with significant dipole strength in the 
XANES of NASB: two due to excitation of the amine core electron (states A and B) and 
one due to the excitation of the nitroso core electron (state C). Given the similarity of the 
XANES spectra, the contributions ESAi and ESA2 to NASB signal (Fig. [HI middle panel) 
are similar to NA. The strongest feature arises due to the double-excitation corresponding 
to states B and C. However, the total signal of NASB at = 806.4 eV is much weaker than 
that of NA, which indicates that these states are uncoupled in the doubly-excited states. 
Analysis of the EGA orbitals (Fig. [7] and [9], middle panels) shows that indeed, the B state is 
strongly localized on the amine group. In NASB the amine and nitroso groups are separated 
by 12.3 A compared to 5.6 A in NA, hence the coupling between the B state (localized on 
amine group) and the G state (localized on the nitroso group) is much weaker than the 
coupling between these states in NA. State A is delocalized, and its coupling with state G 
is significant resulting in the characteristic two-lobe pattern. 

The km spectrum of NADP is very weak (Fig. [HI right panel) with the maximum in- 
tensity approximately 5 times weaker than NA. This is surprising given that the electronic 
stucture of NADP is similar to NASB. Indeed, the equivalent-core orbitals describing the 
singly and doubly core excited states of NADP (Fig. [7] and [HI right panels) are similar 
to those of NASB. The single-orbital approximation would thus predict the 2DXGS signal 
stronger than in NASB due to the shorter distance between the two cores. Analysis of the 
calculated transition dipole moments shows however that the fifejfJ^ejg factor {i.e., the prob- 
ability amplitude of the two-photon excitation into doubly excited state /) is very small due 
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to the cancellation between the transition moments describing excitations via two different 
core shells 



/i/ei/ieig + /i/ea/^eag ~ (13) 

The lowest doubly excited state (marked A+C in Fig. [9l right panel) is described in the 
DFT/ECA approximation by a KS determinant where the + 1 orbital is occupied by the 
nitroso Is electron and A^ + 2 orbital is occupied by the amine electron. The corresponding 



two-photon transition amplitude is given in the EGA by [291] 



where Xno2^ Xnh2 respectively the nitroso and amine N Is orbitals, are 
the valence equivalent- core orbitals describing the promoted core electrons, and <l>(^^\ <l>(^), 
are the Slater determinants made of A^ valence spectator orbitals in four equivalent- 
core molecules. In XANES, the one-electron transition dipole moments (0*^-'^ |/t|xj) are often 
sufficient to qualitatively reproduce the experimental XANES intensities. Factors ($'^-'^|$) 
describe the relaxation among the valence spectator orbitals, which often has a small effect 
on the computed XANES intensities. In case of NADP 2DXCS, however, the calculated 
overlap between the valence orbitals {|<|)(^^) |<|)(2)^ and ^$(^2) |<|)(^)^ have opposite sign, 

($(12)|$(2)^ ^ _^$(12)|^(l)^ (^^5^ 

while the remaining quantities contributing to Eq. (fl^ have the same sign and similar mag- 
nitude. Consequently, the probability amplitude of the corresponding two-photon excitation 
is much smaller than in NA or NASB. This is thus a purely many-body effect induced by 
relaxation among orbitals that are not directly participating in the core transitions. We 
note however that it is difficult to estimate the quality of the EGA in describing many-body 
effects such as the relaxation among the spectator orbitals. More accurate, many-body meth- 
ods that explicitely describe core excitation and, ultimately, comparision with experimental 
data, may be needed to verify the effect of valence relaxation on 2DXCS signals. 
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V. CONCLUSIONS 



The double-quantum-coherence 2DXCS technique proposed here monitors the attosecond 
four-wave-mixing signal in the kjjj phase-matching direction. The contribution of single 
resonances are eliminated and the signal consists only of peaks arising from the coupling 
between the core transitions in the doubly-excited states. Simulations were performed using 
sum-over-states expression derived using the rotating-wave approximation and include the 
pulse envelopes. The differences between km and kj signals were demonstrated on a model 
four-level system with and without coupling between the transitions. Simulations of the 
Nls XANES and kjjj spectra of para-nitroaniline and two-ring hydrocarbons disubstituted 
with an amine and a nitroso group showed that while XANES is virtually invariant to the 
differences in the molecular and electronic structure of these molecules, the double-quantum 
kjii technique is highly sensitive to the separation between the core-shell as well as the 
localization of the corresponding core-excited states. 
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APPENDIX A: 2DXCS SIGNAL OF AN ENSEMBLE OF RANDOMLY- 
ORIENTED MOLECULES 

The first contribution to the km signal [Eq. ([5])] is proportional to 

,,4 3 2 1 _ ai* a-i* a2 ai a4, as a2 ai -ik4r4-ik3r3+ik2r2+ikiri 
H-ge'H'e'fH'feH'eg " ^4 63 6^ IJ'ge'H'e' ff^feH'eg^ i-^^) 

where aj refers to the laboratory-frame Cartesian components of the jth pulse polarization 
vector and the corresponding dipole transition moment, and rj is the laboratory-frame 
position of the core-shell interacting with the jth pulse. 

The transition dipole moments are calculated in the molecular frame and thus must be 
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transformed to the laboratory frame 

Here, /3„ refers to the Cartesian components of the dipole transition moments in the molec- 
ular frame. Iplp^p^p^ is the product of the four directional cosines of the angles between the 
laboratory axes a„ and molecular axes /3„ 

The index / refers to states with two core electrons excited. If e' = e, we have = ri 
and r2 = and 

{gj(fci-fc4)(ri-r2) _ g 

gj(fci-fc3)(ri-r2)^ Yg/ _^ g 

where we used the fact that the 2DXCS signal is measured under the phase matching con- 
dition, fci + A;2 — ^3 — ^4 = 0. 

^-ikir4,-iks,r-i+ik2r2+ikiri _ ^i2-Kq{L/\)cos{9) (A4) 

where 6 is the angle between the vectors ki — ^3 and Vi — r2, q = | fci — fc3|/|A;i|, and L is 
the distance between the two core shells contributing to the signal. Here, 6 depends on the 
orientation of the molecule in the laboratory frame, and q {0 < q < 2) depends on the pulse 
configuration. 



,,4 3 2 1 _ 04 03 ft i27rg(L/A)cos(e)rXXXX \ ( /^rr-s 

The distance between the N and O atoms in NASB, the largest molecule considered, is 
12 A. Assuming g 1, the prefactor is ^^'^l^^'^^W and its contribution to the averaged 

signal is expected to be small. In our calculations we neglected it and used the following 

expression for the XXXX pulse configuration 

,,4 3 21 L r,,/^2 /2 ft ft , ft ft ft ft , ft /3i ft ft] /A«^ 

^J'ge'f'e'ff'fe^^eg " ^5 L^Se'/^e'/Z^/e/^eg + ge' f'e' f fef'eg + /^ge'/^e'//^/e/^eg J l^^J 

Since the third and the fourth pulse have identical polarization, the second term con- 
tributing to Eq. ([5]) is similarly given by Eq. (]A6p . 
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FIGURES 
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FIG. 1: Molecular structure of 4-nitroaniline (NA), 4-nitro-4/-aminestilbene (NASB), and 4-nitro- 
4/-amiiiediphenyl (NADP) 

ii,-^ ,1 ,ik 

k, — 



J.' 



T, T2 T, 



FIG. 2: Pulse sequence in a 2DXCS experiment. Three pulses ki, k2, and ^3 induce a core-hole 
polarization in the molecule, which is probed by the heterodyne field k^. The time intervals ts, t2, 
ti between consecutive pulses serve as control parameters. 
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FIG. 3: Double-side Feynman diagrams representing the two contributions to the kjjj = ki+k2—k3 
signal [Eq. ([6])]. 
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FIG. 4: (Color online) XANES and 2DXCS (imaginary part) of a model four-level system repre- 
senting two uncoupled core transitions. Left panel: The level scheme and XANES. Middle panel: 
The single- and double-resonance contributions and the total kj 2DXCS. Right panel: Two double- 
resonance contributions and the total km 2DXCS. 
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FIG. 5: (Color online) Same as Fig. |3]but for a model four-level system representing two weakly- 
coupled core transitions. 
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FIG. 6: (Color online) Left panel: Experimental ISEELS of nitrobenzene, aniline, and NA from 



Ref. 



45|. Right panel: Experimental ISEELS of NA[45], and present simulated XANES of NA, 



NASB, and NADR 
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FIG. 7: (Color online) Equivalent-core orbitals of NA, NASB, and NADP (as indicated) describing 
the promoted amine or nitroso Is electron in singly-excited states with significant contribution to 
XANES. 
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FIG. 8: (Color online) Two ESA contributions and the total Nls kjjj 2DXCS (imaginary part) of 
NA, NADP, and NASB. 
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FIG. 9: (Color online) Pairs of the equivalent-core orbitals of NA, NASB, and NADP (as indicated) 
describing the promoted amine and nitroso Is electrons in doubly-excited states with significant 
contribution to the corresponding Nls kjjj 2DXCS. 



21 



